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ABSTRACT 


The elastostatic problem for an infinite orthotropic strip con- 
taining a crack is considered. It is assumed that the orthogonal 
axes of material orthotropy may have an arbitrary angular orientation 
with respect to the orthogonal axes of geometric symmetry of the un- 
cracked strip. The crack is located along an axis of orthotropy, 
hence at an arbitrary angle with respect to the sides of the strip. 
The general problem is formulated in terms of a system of singular 
integral equations for arbitrary crack surface tractions. As ex- 
amples Modes I and II stress intensity factors are calculated for 
the strip having an internal or an edge crack with various lengths 
and angular orientations. In most calculations uniform tension or 
uniform bending away from the crack region is used as the external 
load. Limited results are also given for uniform normal or shear 
tractions on the crack surface. 
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1. INTRODUCTION 


Because of the ever increasing use of fiber-reinforced composites 
in a great variety of engineering structures, in recent years the 
problems regarding their structural integrity and failure have been 
studied quite extet!si\ely. In these studies the material is generally 
assumed to be homogeneous and orthotropic if either the structure is 
free from flaws which may be the cause of an eventual failure initia- 
tion, or the structure may have a flaw but its size is large in com- 
parison with the local micros tructural length parameters such as the 
fiber diameter and the distance between the neighboring fibers. Other- 
wise, in failure initiation studies the material has to be treated as 
a nonhomogeneous continuum containing local flaws with certain geome- 
tries. In composites, as well as in wood and certain metallic ma- 
terials, from the viewpoint of structural failure, a distinguishing 
feature of material orthotrOpy is that the material is generally not 
isotropic with respect to its fracture resistance. Furthermore, in 
most cases the planes of orthotropy are generally also the planes of 
weak fracture resistance. Thus, in orthotropic materials regardless 
of the overall geometry and loading conditions, the fracture propa- 
gation would be either along a plane of orthotropy or would have a 
zig-zag path. 

Partly because of the fact that some of the most important 
structural applications of composites have been in sheet form, and 
partly for analytical reasons, the crack problems in orthotropic 
materials have been studied mostly for the cases of plane stress or 
plane strain. In plane problems, if the medium is infinite con- 
taining a line crack or a series of col linear cracks , it was shown 
that the stress intensity factor is identical to that found for an 
isotropic plane with the same crack geometv’y [1-4]. However, it 
was also shown that if the medium is bounded the material orthotropy 
would have an influence on the stress intensity factors, and depend- 
ing on the nature of the orthotropy, the stress intensity factors 
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may be greater or smaller than the corresponding isotropic values [5]. 
In [5] a uniformly loaded orthotropic strip having cracks perpendicular 
to the sides was considered and the plane of the crack was assumed to 
be one of the planes of material orthotropy. This and similar solu- 
tions would be adequate to study the fracture problems in sheet struc- 
tures in which the stress-free boundary is parallel to one of the 
planes of material orthotropy. On the other hand if the stress-free 
boundary of the sheet does not coincide with a plane of orthotropy 
and yet, as expected, if the crack lies on a plane of orthotropy, then 
the solution of the so-called inclined crack would be necessary to 
study the related fracture problem. Such a problem is considered in 
this paper for an infinite strip. The crack is assumed to have an 
arbitrary location and orientation in the strip (Figure 1), the only 
restriction being that the plane of the crack is a plane of material 
orthotropy. The problem is formulated for arbitrary normal and shear 
tractions on the crack surface and the cases of both internal and 
edge cracks are considered. The corresponding internal crack problem 
for an isotropic strip was considered in [6]. 

2. FORMULATION 

The plane elastostatic problem under consideration is described 
in Figure 1 where X] and X 2 refer to the axos of orthotropy and the 
crack is located on the line x-j =0, a<X2<b. The solution of the 
problem is expressed as the sum of two states of stress derived from 
the Airy stress functions F-) (xi .Xg) and F 2 (x,y) where the coordinates 
(xi,X 2 ) and (x,y) are defined in Figure 1. Referred to (x],X2) axes, 
in terms of the stress function Fi the stress components are given by 

( 1 ) ( 1 ) ( 1 ) 

“11 ■ Uf ’ “22 ■ '3x|" ’ “12 ■■ ■ 3x^8Xj ■ ^ ' 

The stress function F^ must satisfy the following differential equa- 
tion [7]: 
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h 3X^3XJ 


IT “ 0 » 


. ,!ll » -!!l2i!66 (3, 

»1 ajj ■ »2 

The elastic constants a^j are defined through the stress-strain rela- 
tions as follows 

Gll = + a ^ 2°22 ’ ^=22 ' ^ 21°11 ^ 22^22 ’ ^^12 ^ 66*^12 


In terms of the engineering constants they are given by 


a^i = 1/Eii , ag 2 = VE22 . ^12 = “^IZ^^IV " ^21 ’ ^66 = ''/‘^12 ' 


By using the Fourier transfonn in the variable X 2 , {2) gives the 
following characteristic equation: 

m‘* - + 0| = 0 . (6) 

Let the roots of the characteristic equation {6) be 


m-| = u)i 


“ -iii2 > ni2 ** ~ '”'^4 


The known constants tiJi and 012 are real if and are complex 

conjugates if . The solution of (2) may then be expressed in 

terms of the following Fourier integrals: 


-u)i s X. 

)e ^ ^ + B(s)e 
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F^(x,^.X2) = 2Tf 


01 , [ S I X, Wp 1 s I X,’ 

At(s)g ^ + B^(s)e 2' 1 


-isxp 

e ‘^ds , x^<0 , 

(8) 


where oi^ and oig are selected in such a way that they have positive 
real parts. Observing that 


Fi(+ 0,X2) = F,(-0,X2) . 8|j-Fi(+<>'X2) =■ • 


(9) 


equations (8) may be written as 


*'i^^r^2^ ■ ir j \ 


A(s)e 


s|x^ 


+ B(s)e 


-oiglslXi 


“iSXr 


'ds , x.j>0 , 


c^A(s) + C 2 B(s) 


0)^ |slx^ 


+ Cc.,A(s) - CiB(s)]e 


a)nls|x,| -1sx« 

'!>e ^ds , x^<0 , (10) 


Where 


'^l " 


o)*! ’to32 


20)r 


2(ii-i 


o),-o)o ' "2 


Cft " 


M '2 


(1), -u 


r“’2 


> c_ = 

3 **C0r 


(11) 


J] '’UJ2 


If we nov/define the discontinuity in the displacement derivatives by 

fgCx) = g-|- c4^ (+0,X2) - U2\-0,X2)] , (12) 

and assume that 

f^(x2) = 0 , fg(x2) = 0 , -«<X2<s , b<X2<“ , (13) 



after soiiie manipulations the unknown functions A and B be obtained 
in terms of f] and f 2 and the stress components may be expressed as 


= 2Tra22(wi^“W^) \ 


( t-X2) f-j ( t)/u^+ii^Xi f2( 


(t-X2)^+w|x^^ 


( i*’'^2 ^ ^1 ^ u)2"tW2Xi f 2 { I') 


( t“X2) 


dt 


JO(v V ) = ] 

27ra22(iJf^i^jr 


t*)l(t-X2)f2(t)-tu^Xi f^(t) 


(t“X2)^+a)^X^ 


^ (t-X2)^+u^X| 


dt , 


a^2\x-,,X2) = 2Trap2(wf-aiS) j 


22^"1 “■2^ Ja 
w^(t“X2)^l ^2^^^ 


{ t-X2 ) f 1 ( t )+u|x^ f 2 ( t ) 
(t-X2)^+wfXj^ 


(t-X2)^+W2^l’' 


dt . 


(14) 


05) 


06 ) 


Referring now to the second solution in which the stress function 
p 2 is expressed in coordinates x,y (see Figure 1), it can be shown 
that the compatibility condition reduces the following differential 
equation: 


a'lFg a'^Fg a'^Fg 3''F2 8“^F2 

Tx^* ‘'^l ax^ay '^2 ax^ay* ^3 ixay^ "'^4 ay‘* ^ 

where 

2Hg 2H^+H3 2Hg _ 

Y-j = - . Y2 = Hg ’ ^3 " ■ "H^ ’ ^4 " 


(17) 


(18) 
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= a^^cos‘*Q + 322510*^6 + ’ 

Hg = a^^s1n'*0 + 322^=°®“® (23^2+^66^®^'^^®^°®^® * 

H3 = agg(cos^0-sin20)^ + 4 ( 3 -| ■j+a22"2a-|2^^^^^®^°^^® ’ 

^4 “ ’’’ ’ 

Hg a [(agg+23i2){cos20-sin'i0) - 2(a^icos20-a22sin20)]sin0coso , 

Hg a [-( 355 + 23 ^ 2 ^'^°®^®“^^'^^®^ " 2(a^isin20-322cos20)]sin0cos0 . 

(19) 

Let the solution of (17) be expressed by the following Fourier 
integral 


p2(x,y) 



p>Co 


4 r,,ys -ixs 

I C^{s)B e ds 


(20) 


Then, substituting from (20) into (17), after some analysis the charac 
teristic equation giving rp..-,r^ is obtained as 

- 173^^ - Y2>'^ ■*■ + 1 = 0 • 

The roots r-j .... .r^ of ( 21 ) are complex and satisfy^ 



( 22 ) 


For the second solution the stress components are found to be 




9y^ 2 it 


4 rf,sy -ixs 

Is^r^C^e/ e ds 
i® 1 


(23) 


(*)fino mAV also note that for 0=0 the roots are real, and if !'i*-->»'4 
are the roots corresponding to the angle 01 , then for ie ang 

0 = 0 r!'ir /2 the roots are ri,..,r 4 . 
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oxDy 27 J f e ds , 

"fti 1 


^ ^ f I s=C„a'’l'-'''e''’‘'ds , 

^(U ^ 


( 24 ) 


(25) 


It will now be assumed that at a given point in the cracked 
orthotropic strip shown in Figure 1 the stress state can be expressed 

nemely equations (14-16) and (23-25), 


<Jtj(^i.X2) ^ J^j.\x^,X2) + cr^^^^(XpX2) > n.J"1.2) , 


(2). 


or 


(1), 


v(x.y) = %gU,y) + o^g(x,y) , („.B) = ^ 


(26) 


(27) 


w?thtr^',v •'" <26) and (27) should be used 

V'flth the following transformations; 


(2) 


”n<'''l.X 2 ) =■ 'lf4x* + "|“yy - SOjUjo'^'' , 

(2) 

yy' 


“latyi-yj) = ("f-"i)Uyy + n,02( 


(2Ba,b) 


and 


Oj^.^(x,y) = + n|crg2 + 2n,n2a,y , 


VTn 


°xy(x,y) = (nf-n|)J/ 2 ^- n^iVgCo/^Ugy) , 
where the direction cosines are given by 


(29a,b) 


111 = COSO , iig = sino . 


(30) 
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3. THE INTEGRAL EQUATION! 


The formulation of the problem given in the previous section 
contains six unknown functions, Cj^(s), (k=l,..,4) and i"j{t), (j=‘1.2). 
Referring to Figure 1, these unknowns can be determined by using the 
following boundary conditions; 

Oyy(x,0) = aj^y(X.O) = <lyy(x,h) = CTj^y{X,h) “ 0 , -«<X<» , (31) 

ffll(0»X2) “ P] ^^2^ * * P2(^2^ ’ a<X2^h (32) 

where the crack surface tractions p-j and P 2 are known functions and 
are assumed to be the only external loads applied to the strip. 
Solutions to other types of loading may be obtained by using the 
standard superposition technique. Substituting from equations (14- 
16, 24, 25, 27, 29) into (31), we obtain the following system of 
algebraic equations expressing C|^(x), k=l,..4, In terms of f] and 
f2: 

4 4 

I C|^(s) = R,(s) . r|,C|,(s) = R 2 (s) , 

4 ri,sh 4 ri,sh 

I C[^(s)e = R 3 (s) , E r^C^(s)e = R^^s) , (33a-d) 

where the functions Rj(s), j=1,..,4, as well aii the solution of the 
algebraic system (33) are given in the Appendix A. 

Substituting now from equations (14, 15, 23-26, 28) into (32) 
and using the appropriate expressions for C[^(s) found in Appendix A, 
the following system of singular integral equations are obtained for 
the functions f-j and fg: 

“^713 I klj(X2.t)]fj(t)dt=p^(X2) . 

a^X2^b , 
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1 = 1,2 , 
(34) 



where 


® * ^2 ** 2 fl 22 (^^'^^ 2 ^ * ( 35 ) 

and the expressions of the kernels k^j, (i»j'=l,2) are given In the 
Appendix B. Referring to the definition of f-| and f£ given by equa- 
tion (12) and the assumptions (13), it is clear that, in addition 
to (34) f'( and f2 must satisfy the following single-valuedness con- 
ditions: 


f.(t)dt « 0 , j=l,2 . (36) 

From the results given in Appendix B, the kernels kij(x2,t), (i,j=l,2) 
appear to be complex valued functions. However, by using the proper- 
ties of the roots w^, (j**l,2) and r|^, (k=l,..,4) of the characteristic 
equations, it can be shown that, as expected, k^j are indeed real 
functions. 

Note that the index of the singular integral equations (34) is 
+1. Therefore, the solution is of the following form: 

f^(t) = q^(t)[(t-a)(b-t)y^ , a<t<b , i»l,2 , (37) 

where the functions g-j and g2 are bounded and continuous in [a,b]. 

It may also be noted that equations (34) give the stress components 
on(0,X2) and o-j2(0,X2) outside as well as inside the region 
(xi = 0, a<X2<b). Therefore, from (34) one may easily obtain the 
stress intensity factors in termc of the unknown functions gi and g2. 
The stress intensity factors are defined by 

k^(a) = lim /2(a-X2)” a^^(0,X2) . 
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(38a-d) 


k2(a) « lim /2(a-"x2) » 

X2->>a 

k-|(b) = lim /2rx"2“b)' a^^(0,X2) » 

Xg-^b 

k2(b) « lim /2(X2-b) <j^2(0»^2^ * 

Xg-^b 


Since the kernels kjjCxg.t), (iJ^'ljE) are bounded in the closed 
interval [a,b], from (37) it follows that the functions 


k^(x2)=^J y. k^j(x2,t)fj(t)dt , (i=l,2) , (O<X2<h/cos0) (39) 

^ a 1 


;b 2 


are also b^.,nJ6U. Thus, defining the fundamental function 


X(z) “ /ri-Fyii^ , (z=X 2 + ix 2 ) (40) 


from (34) and (37) we obtain 





g.(t)dt 


(0=1,2) 


(41) 


Defining now the sectionally holomorphic functions 


*^(z) 



9j(t)dt 

(T-ox^TT) 


(j=l,2) 


(42) 


and observing that <J>i and <>2 are holomorphic outside the cut (a<xg<b, 
X 2 = 0 ), we find 


DjO-jj(x2,0) ” - ‘^j(x2) kj(x2) , (j“l»2, X2’'3 , X2^b) • (43) 

On the other hand, following [8] from (42) it can be shown that 
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94(z) 


( 44 ) 


Where Pj(z) is the principal part of gj/X at infinity. Thus, it is 
seen that 

^ " "X^Cx^V ^ ■*■ * (X 2 <a» X 2 >b) . (45) 

Finally, from (38), (45), and 

X(^ 2 ) = /(X 2 -b)(x 2 -a) = - /(¥^x^y(a^^ , (46) 

we find 

><l (a) = g.|(a)/>^(b-a )/2 , k 2 (a) = -^ 92 (a)/ /(F¥J7^ , 

ki(b) = - -^gi{b)//{b::i72 , l<2(b) = - 92(b)//(b-a)/2 . 

(47a-d) 

4. NUMERICAL SOLUTION AND RESULTS 

The system of singular integral equations (34) is solved numeric- 
ally by fyrst normalizing the interval (a,b) to (- 1 , 1 ) and then using 
the Gauss-Chebyshev integration formulas- [9], The important problem 
in the numerical analysis is the evaluation of the kernels 
To do this a highly accurate and relatively simple technique for the 
calculation of the roots r-j, (i=l,..,4) of the characteristic equation 
(21) was needed. An outline of such a technique may be found in [10]. 
Even though complex algebra had to be used throughout the numerical 
calculations, values of the kernels were, of course* always real. 
First, changing the material constants or the geometry, the isotropic 
results given in [ 6 ] and the results of the symmetric crack geometry 
for the orthotropic strip found in [5] were verified. The numerical 


results are then obtained for the following two basic loading condi- 
tions (see Figure 1): 


- o„,cos20 , p2(x2) * - o^jSinocose , (48a, b) 


which correspond to uniform (membrane) loading +") ~ 

2Xp 

Pl(x2) = Ob^“Ti > 

2x 

PglXg) a COS0-1 )sin0coso (49a, b) 


which correspond to "pure bending." Here ob is the surface stress in 
the strip under bending away from the crack region. Some results are 
also obtained for uniform normal or shear tractions on the crack 
surface in order to explain certain anomalies arising from the in- 
clined crack solution. As an example a boron-epoxy composite sheet 
with the following material constants is considered (see equations 
4 and 5); 


= 24.75 X lOSpsi (170.65 x 10^ N/m^) , 

E22 “ 8 X 106 psi (55.6 X 109 N/m^) , 

G^2 = 0.7 X lOG psi (4.83 X 109 N/m2) , 

v.,2 = 0.1114 . 

For this material the roots m^- or wj, (j=1.2) of the characteristic 
equation (6) turn out to be real. 

The results for the strip containing an internal crack are given 
in Tables 1-4. The stress intensity factors given in the tables are 
defined by equations (38a-d) and are normalized with respect to 


or Ojj/c, c*(b-a)/2. Table 1 shows the results for a syimietric- 
ally located internal crack (i.e., for a» (h/coso) - b) and fbr vari- 
ous values of the angle o. Table 2 shows the results for an excentric- 
ally located internal crack. In this case the crack tip xg^a and the 
crack angle o are fixed (a “ 0.2h/cos8, 0 « t/ 4) and the crack length 
b-a is varied. The stress intensity ratios k-| and kg shown in this 
table are defined in Table 1. The general rule for an excentric crack 
perpendicular to the sides of the strip is that k-|(a) is always greater 
than k(b) if a<h-b. This result is also expected for an inclined crack 
provided the external load is either uniform pressure or uniform shear 
traction on the crack surface. However, in the inclined crack case 
under more general loading conditions this rule may not always be valid. 
For e.xample, from Table 2 it is seen that for b » 0.4h/coso, ki{b)>k^(a). 
Even though this result appears to be somewhat unexpected, it can 
easily be explained by the coupling effect between the shear and nor- 
mal crack surface loadings arising from the inclined crack geometry. 

The stress intensity factors due to only noniial or shear traction on 
the crack surface are shown in Table 3. Note that for the primary 
stress intensity factors (i.e., k-j for normal loading and kg for shear 
loading) the general rule mentioned above remains to be valid. How- 
ever, since the coupling effects (i.e., k^ for shear loading and kg 
for normal loading) can be positive or negative, the type of anomalous 
results observed in Table 2 should not be entirely unexpected. 

In reference [5] it was shown that in an infinite orthotropic 
strip containing cracks perpendicular to the sides the stress state 
in the plane of the crack in general and the stress intensity factors 
at the crack tips in particular are not affected by a 90° rotation of 
the axes of material orthotropy. From the proof given in [53 it can 
be seen that this rather general result will not remain valid for an 
inclined crack. Table 4 shows the result of an example regarding the 
rotation of material axes. In the strip labeled by 30° the stiffer 
material axis E-ji makes 30° with the x-axis, and in that labeled 
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Table 4. Comparison of the stress intensity factors for isotropic 

and orthotropic strips with a symmetrically located internal 
crack. Tension: » bending: oyu*= ok(l-2y/h) , 

(b-a)/{h/cos0) = O.b, c=(b-a)/2» a = (h/cos0)-b 



<X> 

II 

O 

0 = ir/6 

Tension 

Tension 

Bending 



h'%-^ 

ki/a[j/c 


Isotropic 

1.303 

1.080 

0.504 

0.248 

0.137 

Ortho. (30“) 

1.226 

1.420 

0.553 

0.288 

0.141 

Ortho. (120®) 

1.226 

1.172 

0.518 

0.258 

0.138 


by 120® Eti axis makes 120“ with the x-axis, i.e., in the latter case 
the material has been rotated by 90“ (see Figure 1). The isotropic 
results are also given in the table. The table shows that in the in- 
clined crack problem not only the material orthotropy but also the 
orientation of the axes of orthotropy may have a significant effect 
on the stress intensity factors. 

In the case of an edge crack, i.e , for a = 0, b<h/cos0, the in- 
tegral equations (34) remain unchanged. H'',w«;ver> the unknown functions 
fl (t) and f2(t) are bounded at t = 0 and thti conditions (36) are no 
longer valid. In this case the integral equations can be solved nu- 
merically by first normalizing the interval (0,b) to (-1,1) through 
the change in variables 

t = |(r+l) , X2=|(s+1) , -l<(s,r)<l , (50) 

and then using again a Gauss-Chebyshev integration formula. A conven- 
ient technique in this problem is defining the unknown functions by 

f^(t) = G.(r)//Tr^ , i=l,2 (51) 
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and using the collocation points sj obtained from U^„]{sj) = 0, 
(j=li..,n-l) and the condition G-j(-l) = 0 (to account for boundedness 
of f-j(t) at t = 0) to calculate 6-j(ri^), (k«l,..,n) Tn(r|<) = 0, where 
Tji and are Chebyshev polynomials. Table 5 shows the calculated 
results for the edge crack. In this problem too the external load 
is either a uniform tension or a uniform bending applied to the strip 
away from the crack region. 

ACKNOWLEDGEMENTS 

This work was partially supported by the National Science 
Foundation under the Grant ENG77-19127, by NASA under the Grant 
NGR39-007-011 and by the Technical University of Istanbul. The 
numerical calculations were carried out at the Computing Center 
of the Technical University of Istanbul . 


- 17 - 




OP POOR QUALITY 

- 18 - 







REFERENCES 


], Ang, D. D. and Williams, M. L., "Combined stresses in an ortho- 
tropic plate having a finite crack," Journal of Applied Mechanics 
Vol. 28, Trans. ASME, 1961, pp. 372-37F; 

2. Savin, 6. N., Stress Concentration Around Holes . Pergainon Press, 
New York, 1961. 

3. Sih, G. C., Paris, P, C. and Irwin, G. R. , "On cracks in recti- 
linearly anisotropic bodies," Int. J. Fracture Mechanics , Vol. 1, 
1965, pp. 189-203. 

4. Krenk, S., "Stress distribution in an infinite anisotropic plate 

with col linear cracks," Int. J. Solids and Structures, Vol. 11, 
1975, pp. 449-460. “ 

5. Delale, F. and Erdogan, F., "The problem of an internal and edge 
cracks in an orthotropic strip," Journal of Applied Mechanics, 
Vol. 44, Trans. ASME, 1977, pp. 237-242. 

6. Krenk, S., "On the elastic strip with an internal crack," Int . 

J. Solids and Structures . Vol. 11, 1975, pp. 693-708. 

7. Lekhnitskii, S. G. . Theory of Elasticity of an Anisotropic 
Elastic Body , Holden-Day*. Inc. 1963.* 

8. Muskhelishvili , N. I., Singular Integral Equations , P. Noordhoff, 
Groningen, The Netherlands, 1953. 

9. Erdogan, F. and Gupta, G. D., "On the numerical solution of singu 
lar integral equations," Quarterly of Applied Mathematics , Vol. 
30, 1972, pp. 525-534. 

10. Delale. F. and Erdogan, F., "Transverse shear effect in a circum- 
ferentially cracked cylindrical shell," NASA Technical Report, 
Grant No. N6R309-007-011, July 1977. 


glUGlNAL PAGE IS 
OF POOB QUiU^rrY 


- 19 - 


APPENDIX A 


(33): 


Expressions of the functions Rj(s) and the solution of equations 


n$ 

fc\ T 1 / C. 


R^{s) a - J^)lj + (o)|-f^fwi)l 2 + (nfoij -n|wi)Jj 

+ (n^wg - + En^ngCw-iKj -WgKj-w^Lj +0)21^]} 

+ ( n^-Ug) [^1 Kj ~ WgKg - u)i L| + WgLg]} j 

n ^ n ^ 

R3(s) = A^{(nftu^ ■*■ ■*■ 

+ (n^wg - n|oj^)J^ + 2n^ng[(o^K^ -UgKg-w^^ +{DgL^]} 
^4(5) = ^riTn 2 [“'(jjj-+u)^)I| + (t^+Wg)I| - (u)^+{])|)J| 

+ {u)2+ui|)J^] + (n^^“n2)[w^Kf-uigK^-(D^L|+w2L|]} ; 

^’1 ” 2TOg2(w|'-a)|)s'‘ ’ 

Ij{s) = f E^(5,t)f,(t)dt , 

jk(5) , F5(s,t)f2(t)dt , 

^ Ja 

Kj!s)=fF5(s,t)f,(t)dt , 

L§(s) = |VE 5 (s.t)f 2 (t)dt , (j,k)“( 1 . 2 ) ; 
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,XJC0S Cjt . c^sInHcjt 


"IsltrinX. 

Ej(s,t) » ire *^[n 

+ 1 jlp CjCOS CjSt - iriiXjSin c^st] , 
r2, .V -ls!x^(h-n,t) 

^j{s,t) = ue {'■niXjCos[cjs(t-n^h+n^wjh)] 

+ CjSin[|s( Cj{t-n^h+n^wjh)] 

+ 1cj jly cos [scj(t-n^h+wjn^h)] 

+ in^XjSln [scj(t-n^h+wjn^h)]} , j«l,2 ; (AlO.All) 

1 “Islx^n^t Cj 

Fj{s,t) = e [- cos CjSt + n-jbjSin Cj|s|t 

J 

+ i -[Tj- n^bjCos CjSt + 1 ^ sin CjSt] , 

2, ”ls|X.(h-nit) Cj 

Fj(s,t) =ne ^ cos[scj(t-n^h+wjn^h)] 

+ n^bjSin[ js|Cj(t-nih+aj^n^h)] 

+ i -||-j-nibjCOs[sCj(t-n^h-HDjnih)] 

c,- 

- i ojt sin[scj(t-n|h+ujn|h)] , j-1,2 ; (A12.A13) 

Xj " Wj/(n|aij+n|) , bj = V(n|(oj+n|) , Cj = -n2/(n|w?+n|) , 


J"1 >2 j 


(A14-AVj) 


Solution of equations (33); 
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r^sh Tgsh r^ish r^sh 


A(s) « ' -e )(e ^ - e ) 

ToSh PoSh r-|Sh r^sh 

- (ri“P3)(rpr2){e ‘- -e ^ )(e ' -e^ ) 

PoSh p^sh p,sh p-sh 

+ (ri-r* 2 )(p^'P 4 Me ^ -e ^ )(e -e ^ ) ; 

(pp+p,)sh (po+»'d)sh 


+ r2(r^-r3)e 


(f3+l-4)sh 


(p„+p-)sh (p«+p,)sh 

»ni2(s) = {P2“»"3)e ^ ^ - (P2“p4)e ^ ^ - (p4~P3)' 


V sh sh r* s 

0113 ( 5 ) » 1'2('"4‘''3^® ^ ^ ^ 

PnSh PgSh p^sh 

11114(5) = (i’3-i’4)e + (p4“P2)e + (>'2'‘'3^® ’ 


m 


21 


(s) = P4(Pi"P3)e 


(Pl+p3)sh 

{p3+r^)sh 




(Pl+P4)sh 


+ Pl(P3-P4)e 

(pi+P3)sh 

0122 ( 5 ) = (i'i-i'3)e + (^T*'4)® 

(P3+P4)sh 

+ (P4-P3)e 


(Pl+p^)sh 


PlSh 


PsSh 


0123(5) = ri(P3-P4)e - p3(Pi-p4)e + l'4(''‘r^3^® 


P 4 SI 
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(A17) 


(A18) 


(P3+P4)sh 



/ . . ^nsh r>,sh 

ni 24 (s) = (r 4 -r 3 )e + (r^-r^)e ^ ^ 


1 ^ 31 ( 5 ) « (r 4 -r 2 )e^'" 2 ’^‘'^ 


)sh 


+ (r^-r^)e 


(»’l+r 4 )sh 


ni 33 {s) = r^(r^-r 2 )e - r 2 (r^-r 4 )e - r 4 (r^-f 2 )e ^ 

n^34(s) = (r2-r4)e - ir■^-r^)e + (V'^-r2)e 

. (r^+rglsh (r,+r 3 )sh 

ni4^(s) = r3(r^-r2)e ' ^ 


+ r^(r 2 -r 3 )e 


(r 2 +r 3 )sh 


m 


)sh 


'42 


, . (r,+rp)sh (ri+r, 

(s) = - (r^'>r2)e ^ + (^^.^3)0 ^ 

+ v)^ 3 -r 2 )e , 


/X , V r,sh 

11143(5) = r^(r2 -r3 )e - ■*■ '"3(^r*’2^® 


"’ 44(0 « (»' 3 -f' 2 ^® V ■‘* 


r/jsh r,sh 

- (r,-C 2 )e 


(A 19 -A 34 ) 
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APPENDIX B 


ORIOINAI' PA<iK IS 
OF puui: uI’ai.h y 


Expressions of the kernels 


k,,(,,.t) . d, |■[S,J(x2.t.s) * G,j(«2.t.-snds . (1.4) ■ (1.2) 1 

(Bl) 


'0 




• 102 X 2 $ 


6u<’‘2'‘'‘> • HTsT t)E{‘'’2Wi*'’4M’'V''lM 

- 1l5j2F‘th5U,F{-‘.5U2Fi] • 


- 102 X 2 $ 


G,2(.2.t.F) ■ Hjrr [• '^'’5E|*u2h5Ej-u,h5E{*...2hjE| 

♦ u)Jh^F{+u^h2F^-^^|h2F{+w^h^F| 

P,E{*V2EJ- ;E{*VjE|*V5..^F{ 

. VjuijFj-fVju^Fj-VjUjFlJ . 


^21 


- 102 X 2 $ 


(B2.B3) 


G22(x2.t.F) ■ ‘ a( ?T [--i'' 5^'*“2’'5^*‘1''6^1*^2''64 
♦ -;v,F{*-2-.2fH"3''f*“2''4''2] '' 


Where the functlohs Ej(x.t) end Fj(s.t). (4.E-1.2) ere given by eque- 
tions (A10-A13K ^($) 1$ Q^vcn by (A18), and 
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4 


h^(Xj,.s) ■ 

I <X|^(x2>s) ’ 


h2(x2,s) ■ 

4 


h3(x2,s) - 

4 

1 ^k^*l‘"k3'*’^*A4^ * 


h4(x2,s) - 

4 

j “k^*A3^^*4'"k4^ • 


hg(x2,s) - 

4 

1 “k[2n^n2in,^^-i(n^^-n|)m,^23 

• 

hg(x2.s) ■ 

4 

f \C2n^"A3-^("f-'’2K43 

; (B8-B13) 

\{x2.s) • 

, r.n,x-s 

(n{rj-ri^-21n^n2rk)e ^ 

. {i-l...,4) ; (E’.4) 

v^(x2,s) • 

4 

f Sj(aiuiji+1a2mj2) . 


V2(x2,s) - 

4 

f 6j(a3"’jl+<Vj2^ ’ 


V3(x2,s) - 

4 

1 Sj(*i"’j3‘^^32mj4) , 


V4(x2.s) - 

4 

[ 6j(a3m.3Ha4mj4) . 


V5(x2,s) ■ 

4 

1 0j[2n^n2.Tij^-1(n?-n^)mj2] 

1 

4 

[ 6j[2n^n2mj3-1{n{-n^)mj4l 

, (B15-B20) 

Vg(x2,s) « 

» 

6j(x2.s) 

• [n^n2rj+nin2^1(n{-n^)rj]e ^ ^ ^ , (j»l,..,4) (B21) 
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(B22-B25) 


*1 ’ "i“i * ’ *2 " ~ ^ "i"2‘^i ’ 

nl n ^02 

«3 " '"f‘^‘2 * ^ ’ ^4 " ‘ ’ "l"2‘**2 ’ 

and the functions %j(s). (k,J-l,..,4) are given by equations (A19-A34). 




